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Abstract In this paper, we constrain the Cardassian expansion models from the lat- 
est observations including the updated Gamma-ray bursts (GRBs), which calibrated 
cosmology-independently from the Union2 compilation of type la supernovae (SNe la). 
By combining the GRB data to the joint observations with the Union2 SNe la set, 
along with the Cosmic Microwave Background radiation observation from the seven-year 
Wilkinson Microwave Anisotropy Probe result, the baryonic acoustic oscillation obser- 
vation from the spectroscopic Sloan Digital Sky Survey Data Release galaxy sample, 
we find significant constraints on model parameters of the original Cardassian model 
Quo = 0.282t°;°^, n = 0.03±g;g|; and n = -0.16^1, P = 0.76tg;^ of the modified 
polytropic Cardassian model, which are consistent with the ACDM model in 1-er confi- 
dence region. From the reconstruction of the deceleration parameter q(z) in Cardassian 
models, we obtain the transition redshift zt = 0.73 ± 0.04 for the original Cardassian 
model, and zx = 0.68 ± 0.04 for the modified polytropic Cardassian model. 

Key words: Gamma rays : bursts — Cosmology : cosmological parameters 
1 INTRODUCTION 

Recent years, the cosmological observations from type la supernovae (SNe la; Riess et al.1998; 
Perlmutter et al. 1999), cosmic microwave background radiation (CMB; Spergel et al. 2003), and large- 
scale structures (LSS; Tegmark et al. 2004; Eisenstein et al. 2005) have been used to explore cosmology 
extensively, which support that the present expansion of our universe is accelerating. In order to explain 
the accelerating expansion of the universe, many cosmological models have been proposed. The first 
categories are proposed by introducing an energy component called dark energy with negative pressure 
in the universe, which dominates the universe to drive the acceleration of expansion at recent times. 
Many candidates of dark energy have been taken into account, such as the cosmological constant with 
equation of state w = — 1 (Carroll et al. 1992), the scalar field models with dynamical equation of state, 
e.g., quintessence (Ratra & Peebles 1988; Caldwell et al. 1998), phantom (Caldwell 2002), k-essence 
(Armendariz-Picon et al. 2001), tachyon (Padmanabhan 2002, Sen 2005), quintom (Feng et al. 2005; 
Guo et al. 2005; Liang et al. 2009); as well as the Chaplygin gas (Kamenshchik et al. 2001) and the 
generalized Chaplygin gas model (GCG, Bento et al. 2002), the holographic dark energy (Cohen 1999; 
Li 2004), the agegraphic dark energy (Cai 2008; Wei & Cai 2008a, 2008b), the Ricci dark energy (Gao 
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et al. 2009) and so on. On the other hand, many alternative models in which gravity is modified to drive 
the universe acceleration have been proposed, e.g., the f(R) theory in which the non-linear gravity 
Lagrangian (L ~ R + f(R), where R is the scalar curvature) has been taken into account (Capozziello 
& Fang 2002); the braneworld models such as the Dvali-Gabadadze-Porrati (DGP) model which con- 
sider that our observable universe might be a surface or a brane embedded in a higher dimensional 
bulk spacetime (Dvali et al. 2000); as well as the Cardassian expansion model in which the Friedmann 
equation is modified (Freese and Lewis 2002; Wang et al. 2003). 

In 2002, Freese and Lewis (Freese & Lewis 2002) proposed the Cardassian expansion model as a 
possible explanation for the acceleration by modifying the Friedmann equation without introducing the 
dark energy. The modified Friedmann equation for the original Cardassian model is 

The Cardassian term which is proportional to p n may show that our observable universe as a 3 + 1 
dimensional brane is embedded in extra dimensions. The first term on the right side of the equation 
dominates initially, so the equation becomes to the usual Friedmann equation in the early universe. 
Then the two terms become equal at redshift z = z ca rd ~ 0(1) (Freese & Lewis 2002), and thereafter 
the Cardassian term begins to dominate the universe, n is assumed to satisfy n < 2/3 to give rise 
to a positive acceleration of the universe. If n = 0, the Cardassian term becomes the cosmological 
constant. If B = 0, the equation becomes the usual FRW equation without the cosmological constant. 
Furthermore, the modified polytropic Cardassian model can be obtained by introducing an additional 
parameter f3 into the original Cardassian model (Wang et al. 2003), and the corresponding modified 
Friedmann equation is 

H 2 = 8nG^ + Cp np ) i/p^ (2) 

When (3=1, the modified polytropic Cardassian model reduces to the original model. 

Gamma-ray bursts (GRBs) are likely to occur in high-redshift range beyond the SNe la redshift 
limit. Up to now, the farthest GRB detected is GRB 090423 at z = 8.2 (Tanvir et al. 2009; Salvaterra 
et al. 2009). Recently, several empirical GRB luminosity relations have been proposed as distance indi- 
cators (Amati 2002; Norris et al. 2000; Fenimore & Ramirez-Ruiz 2000; Reichart et al. 2001; Schaefer 
2003a; Yonetoku et al. 2004; Ghirlanda et al. 2004a; Liang & Zhang 2005; Firmani et al. 2006a; Yu et 
al. 2009). Therefore, GRBs could be regarded as the standard candles to be a complementary cosmo- 
logical probe to the universe at high redshift (Schaefer 2003b; Takahashi et al. 2003; Bloom et al. 2003, 
Dai et al. 2004; Ghirlanda et al. 2004b; Friedman & Bloom 2005; Firmani et al. 2005, 2006b; Liang & 
Zhang 2005; Di Girolamo et al. 2005; Bertolami & Silva 2006; Ghirlanda et al. 2006; Schaefer 2007; 
Wright 2007; Wang et al. 2007; Amati et al. 2008; Basilakos & Perivolaropoulos 2008; Cuesta et al. 
2008a, 2008b; Daly et al. 2008; Qi et al. 2008a, 2008b; Vitagliano et al. 2010). Due to the lack of the 
low-redshift sample, these luminosity relations have been usually calibrated by assuming a particular 
cosmological model (e.g., the ACDM model with particular model parameters according to the con- 
cordance cosmology). Therefore most of the calibration of GRBs are always cosmology-dependent to 
lead the circularity problem in cosmological research. Many of works treat the circularity problem with 
statistical approaches such as simultaneous fitting of the parameters in the calibration curves and the 
cosmology (Li et al. 2008; Wang 2008; Samushia, & Ratra 2010; Xu 2010; Graziani 2010). However, it 
is noted that an input cosmological model is still required in doing the simultaneous fitting. 

In our previous paper (Liang et al. 2008), we presented a new method to calibrate GRB luminosity 
relations in a completely cosmology-independent way. The luminosity distance of GRBs in the redshift 
range of SNe la can be obtained by interpolating directly from the SNe la Hubble diagram, and GRB 
data at high redshift can be obtained by utilizing the calibrated relations (Liang et al. 2008). Similar 
to the interpolation method, the luminosity distance of GRBs could be obtained by other mathematical 
approach, such as the empirical formula fitting (Kodama et al. 2008), the non-parametric reconstruction 
method (Liang & Zhang 2008), the local regression (Cardone et al. 2009), and the cosmographic fitting 
(Gao et al. 2010; Capozziello & Izzo 2010). Following the GRB calibration method directly from SNe 
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la, the derived cosmology-independent GRB data at high redshift can be used to constrain cosmological 
models by using the standard Hubble diagram method (Liang & Zhang 2008; Capozziello & Izzo 2008; 
Izzo et al. 2009, Izzo & Capozziello 2009; Wei & Zhang 2009, Wei 2009; Wang et al. 2009a, 2009b; 
Qi et al. 2009; Wang & Liang 2010; Liang et al. 2010; Wei 2010, Freitasa et al. 2010; Liang et al. 
2011; Liang & Zhu 2011). Capozziello & Izzo (2008) first used the GRB relations calibrated with 
the so-called Liang method to derive the cosmography parameters at high redshift. Liang et al. (2010) 
combined the GRB data with the joint data to constraint the cosmological parameters and reconstructed 
the acceleration history of the universe. 

Here we consider the Cardassian model viewed as purely phenomenological modifications of the 
Friedmann equation to drive the universe acceleration and focus on the latest cosmological constraints 
including GRBs. Up to now, the Cardassian model have been constrained from many observational data, 
such as the angular size of the compact radio sources (Zhu & Fujimoto 2002), SNe la (Wang et al. 2003; 
Zhu & Fujimoto 2003; Szydlowski & Czaja 2004; Godlowski et al. 2004; Frith 2004; Bento et al. 2005), 
the x-ray gas mass fraction of clusters (Zhu & Fujimoto 2004; Zhu et al. 2004), CMB (Sen & Sen 2003; 
Savage et al. 2005), the large scale structure (Multamaki et al. 2003; Amarzguioui et al. 2005; Fay & 
Amarzguioui 2006), the gravitational lensing (Alcaniz et al. 2005), the baryonic acoustic oscillation 
(BAO) (Wang et al. 2007), the Hubble parameter versus redshift data (Yi & Zhang 2007), as well as the 
different combined data (Bento et al. 2006; Davis et al. 2007; Wang 2007; Wang & Wu 2009; Feng & 
Li 2010). Also, constraints from GRBs with the joint analysis on the Cardassian model can be obtained 
in (Wang et al. 2007; Cuesta et al. 2008a; Wang et al. 2009a; Wang & Liang 2010). Very recently, the 
Union2 compilation of SNe la data set which consists of 557 SNe la has been released (Amanullah et 
al. 2010), whereas the seven-year data of Wilkinson Microwave Anisotropy Probe (WMAP7) has also 
been released (Komatsu et al. 2010). In this paper, with the updated GRB data calibrated directly from 
the Union2 set, we constrain the Cardassian model and the modified polytropic Cardassian model from 
the latest observations by combining the GRB data to the joint observations with the Union2 set, along 
with the CMB observation from the WMAP7 result, the BAO observation from the spectroscopic Sloan 
Digital Sky Survey (SDSS) Data Release galaxy sample (Eisenstein et al. 2005). We also reconstruct 
the deceleration parameter q(z) in Cardassian expansion models and obtain the transition redshift zt. 
We find that tighter and more stringent constraints can be given out with the combined data including 
GRBs in this work. 

This paper is organized as follows. In section 2, we introduce the analysis for the observation data. 
In section 3, we present constraint results on Cardassian models from the joint observations including 
GRBs, as well as SNe la, CMB, and BAO. Conclusions and discussions are given in section 4. 



2 OBSERVATIONAL DATA ANALYSIS 

In our previous paper (Liang et al. 2008; Liang et al. 2010), we used the 192 SNe la compiled by Davis 
et al. (2007) and the 397 SNe la set (Hicken et al. 2009) in the interpolation procedure to calibrate GRB 
luminosity relations from the 69 GRBs compiled by Schaefer (2007). A larger number of SNe la sample 
could bring more accurate result in the interpolation procedure. Very recently, the Union2 compilation 
(Amanullah et al. 2010) of 557 SNe la data set has been released by the Supernova Cosmology Project 
Collaboration (SCP). In this paper, we use the Union2 set to calibrate GRB luminosity relations with 
the GRB sample at z < 1.4 by using the linear interpolation method, and we update the distance moduli 
of the GRBs at z > 1.4 obtained by utilizing the new calibrated relations. For more details for the 
calculation, see Liang et al. (2008); Liang et al.(2010). We plot the Hubble diagram of Union2 SNe la 
and the GRBs obtained using the interpolation methods in figure 1 . The distance moduli of the 27 GRBs 
at z < 1.4 are obtained by using the linear interpolation method directly from the Union2 SNe set; the 
42 GRB data at z > 1.4 are obtained by utilizing the five relations calibrated with the sample at 2 < 1.4. 

The position of the first acoustic peak in the power spectrum of CMB favors a spatially flat Universe, 
therefore we assume a flat universe prior throughout this work. Constraints from SNe la and GRB data 
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Fig. 1 Hubble Diagram of 557 SNe la (red dots) and the 69 GRBs (circles) obtained using the 
interpolation method. The 27 GRBs at z < 1.4 are obtained by linear interpolating from SNe 
la data (black circles), and the 42 GRBs at z > 1.4 (blue circles) are obtained with the five 
relations calibrated with the sample at z < 1.4. The curve is the theoretical distance modulus 
in the concordance model (w = — 1, Hmo = 0.27), and the vertical dotted line represents 
z = 1.4. 



can be obtained by fitting the distance moduli A distance modulus can be calculated as 

/i = 5 log + 25 = 5 log 10 Dl — Mo, (3) 

where /l*o = 5 log 10 h+42.38, h = Ho j (lOOkm/s/Mpc), Hq is the Hubble constant. For a fiat universe, 
the luminosity distance Dl can be calculated by 



D L = H d L = (l + z) ——, (4) 
Jo E(z') 

where E{z) = H j Ho, which determined by the choice of the specific cosmological model. The \ 2 
value of the observed distance moduli can be calculated by 

^ = 1^ ~2 . (5) 

i=l A*>* 

where /Uobs(-Zj) are the observed distance modulus for the SNe la and/or GRBs at redshift z\ with its 
error <7 W ; fi(zi) are the theoretical value of distance modulus from cosmological models. Following the 
effective approach (Nesseris & Perivolaropoulos 2005), we marginalize the nuisance parameter no by 
minimizing 



Xl = A - B 2 /C, (6) 
where A = J2 \poba(zi) - 5 log 10 Drf/a 2 ^, B = £ \fJ, ohs (zi) - 5 log 10 Dl]/^, and C = £ 1/cr^. 
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For the CMB observation, the shift parameters R provide an efficient summary of CMB data to 
constrain cosmological models. For a flat universe, the shift parameter can be expressed as (Bond et al. 
1997) 

R = / ~EFT~\~ ( ? > 

Jo E(z) 

here z rcc is the redshift of recombination. From the WMAP7 result, the shift parameter is constrained 
to be R = 1.725 ± 0.018 and z Iec = 1091.3 (Komatsu et al. 2010). The \ 2 value of the shift parameter 
can be calculated by 

2 _ (i?- 1.725) 2 

XCMB ~ q q 18 2 ■ (°> 

For the BAO observation, we use the distance parameter A which can be expressed as for a flat 
universe (Eisenstein et al. 2005) 



E(z BAO ) l / 3[ z BAO J E(z) 1 K> 

where zbao = 0.35. From the SDSS spectroscopic sample of luminous red galaxy, the distance pa- 
rameter is measured to be A = 0.469(n s /0.98)~ ' 35 ± 0.017 (Eisenstein et al. 2005) , with the scalar 
spectral index n s = 0.963 from the WMAP7 data (Komatsu et al. 2010). The \ 2 value of the distance 
parameter can be calculated by 

2 _ (A - 0.4666) 2 
Xbao ~ 0.017 2 ' 

3 CONSTRAINTS FROM COMBINING GRBS, SNE IA, CMB, AND BAO 

In order to combine GRB data into the joint observational data analysis to constrain cosmological mod- 
els, we follow the simple way to avoid any correlation between the SNe la data and the GRB data (Liang 
et al. 2010). The 40 SNe points used in the interpolating procedure are excluded from the Union2 SNe 
la sample used to the joint constrains. The best fit values for model parameters can be determined by 
minimizing 

X = X A1 {5i7SNc+42GRBs} + XCMB + XBAO ■ (H) 

From the modified Friedmann equation of the original Cardassian model, if only considering the 
matter term without considering the radiation for simplification, using pu = Pmo(1 + -z) 3 = S7mo/0c(1 + 
z) 3 , where the present critical density of the universe p c = 3Hq/8ttG, we obtain 

f x (z) = £*- = (l + z) 3 ", (12) 
Pxo 

where subscript 'x' refers to any component providing additional term in the Friedmann equation. The 
corresponding E(z) of the Cardassian model is 

E(z) = [f2 M0 (l + zf + (1 - fi M o)(l + z) 3 "] 1/2 . (13) 
For the modified polytropic Cardassian model, we obtain 

The corresponding E(z) of the modified polytropic Cardassian model is 

E(z) = [Q M0 (1 + z?[l + J0^ ¥ ] 1, "\ 1/2 d5) 
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Original Cardassian Model 


Modified polytropic Cardassian Model 
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Table 1 The best-fit value of parameters Omoi n an d ft for the original Cardassian model and 
the modified polytropic Cardassian model with la uncertainties, as well as Xmin> Xmin/d°f> 
with SNe+GRBs+CMB+BAO, SNe+CMB+BAO, and GRBs+CMB+BAO. 




Fig. 2 The joint confidence regions in the £1mo — n plane for the original Cardassian 
model in a flat universe. The contours correspond to 1-er and 2-a confidence regions. 
The black solid lines, red dashed lines, the blue dash-dotted lines, represent the results of 
SNe+GRBs+CMB+BAO, SNe+CMB+BAO and GRBs+CMB+BAO, respectively. The black 
plus, red point, and blue star correspond the best-fit values of SNe+GRBs+CMB+BAO, 
SNe+CMB+BAO and GRBs+CMB+BAO, respectively. 



The joint confidence regions in fluo — n plane with the combined observational data for the original 
Cardassian expansion model are showed in figure 2. With SNe la + GRBs + CMB + BAO, the best-fit 
values at 1-er confidence level are Qmo = 0.282^ooi4, n = 0.03^o'o5- F° r comparison, fitting results 
from the joint data without GRBs are also given in figure 2. The best-fit values with SNe la + CMB + 
BAO are Sl M n = 0.270^;^, n = O.CJO+^5. While with GRBs + CMB + BAO without SNe la, the 
best-fit values are Hmo = 0-290^o'o46' n = O.ll^Q'25. We present the best-fit value of f^Mo* n with \-a 
uncertainties, and xLi a > Xmin/^of for the original Cardassian model in Table 1. 
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Fig. 3 The joint confidence regions in the n — plane for the modified polytropic Cardassian 
model in a flat universe. The contours correspond to 1-er and 2-er confidence regions. 
The black solid lines, red dashed lines, the blue dash-dotted lines, represent the results of 
SNe+GRBs+CMB+BAO, SNe+CMB+BAO and GRB s+CMB+B AO, respectively. The black 
plus, red point, and blue star correspond the best-fit values of SNe+GRBs+CMB+BAO, 
SNe+CMB+BAO and GRBs+CMB+BAO, respectively. 



For the modified polytropic Cardassian model, we find that the best-fit values at 1 -a confidence level 
with SNe la + GRBs + CMB + BAO are O M0 = 0.285^;^, n = -0.16±£^, (3 = 0.76+^f. Figure 
3 shows the joint confidence regions with the combined observational data for the modified polytropic 
Cardassian model in the n — f3 plane, while fixing J7mo with the best-fit values. With SNe la + CMB + 
BAO, the best-fit values are Q M0 = 0.271+^u, n = -0.22^°;^, /? = while with GRBs 

+ CMB + BAO, the best-fit values are fl M o = 0.285±£^, n = -0.06^;^, /3 = 0.81^;^. We also 
present the best-fit value of Omo. n, and f3 with \-a uncertainties, and Xmin< Xmin/ ^of for the modified 
polytropic Cardassian model in Table 1 . 

From Fig. 2-3 and Table 1, we can find that GRBs can also give strong constrains when combined 
to CMB and BAO data without SNe la. By comparing to the joint constraints with GRBs and without 
GRBs, we can see that the contribution of GRBs to the joint cosmological constraints is a slight shift 
which adding the best-fit value of QmO) an d significantly narrowing the parameters confidence ranges 
of the modified polytropic Cardassian model. We also find that the ACDM model (n = 0, /3 = 1 ) 
is consistent with all the joint data in l-a confidence region, and combining these observational data 
can tighten the model parameters significantly comparing to results from former works (Wang 2007; 
Wang & Wu 2010). We also investigate the deceleration parameter for Cardassian expansion models. 
The deceleration parameter q(z) can be calculated by 



x , . i dE(z) 
q = -l + (l + z)E(z)- 1 —y- 



(16) 
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Fig. 4 The evolution of the deceleration parameter q(z) from fitting results in the original 
Cardassian model. The solid line is drawn by using the best fit parameters. The shaded region 
shows the l-a uncertainties. 

In figure 4, we show the evolution of q(z) for the original Cardassian expansion model. We obtain 
g = —0.55 ± 0.054, and the transition redshift is zx = 0.73 ± 0.04 at the lcr confidence level, which 
is more stringent and comparable with the former result (zt = 0.70 ± 0.05) by Wang (2007), but is 
slightly later than the former result (zt = 0.55 ± 0.05) by Wang & Wu (2010) . We show the evolution 
of q{z) for the polytropic Cardassian expansion model in figure 5, and we find the transition redshift 
z T = 0.68 ± 0.04 and q = -0.57 ± 0.07. 

4 CONCLUSIONS AND DISCUSSIONS 

In this paper, by using the Union2 set of 557 SNe la, we calibrate GRB data in a completely cosmology- 
independent way. When combine the GRB data with the Union2 set, we avoid any correlation between 
the SNe la data and the GRB data (Liang et al. 2010). From GRB data to the joint observations with 
the Union2 set, along with the CMB from WMAP7 and the BAO observation from the SSDS Data 
Release galaxy sample, we find significant constraints on model parameters of the original Cardassian 
model n M o = 0.282^;°^, n = 0.03^;^; and n = -0.16±£^, f3 = 0.76^;^ of the modified 
polytropic Cardassian model, which are consistent with the ACDM model in l-a confidence region. 
From the reconstruction of the deceleration parameter q(z) in Cardassian expansion models, we obtain 
the transition redshift zt = 0.73 ± 0.04 for the original Cardassian model, and zt = 0.68 ± 0.04 for the 
modified polytropic Cardassian model, which are more stringent comparing to the former results (Wang 
2007; Wang & Wu 2010). It is found that GRBs can give strong constrains when combined to CMB and 
BAO data without SNe data, and we can see that the contribution of GRBs to the joint cosmological 
constraints by comparing to the joint constraints with GRBs and without GRBs. Hereafter, along with 
more and more observed data, GRBs could be used as an optional choice to set tighter constraints on 
the Cardassian model and even other cosmological models. 
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Fig. 5 The evolution of the deceleration parameter q(z) from fitting results in the poly tropic 
Cardassian model. The solid line is drawn by using the best fit parameters. The shaded region 
shows the \-a uncertainties. 



Recently, some works point out that there are observational selection bias in GRB relations (Butler 
et al. 2007; Shahmoradi & Nemiro 2009) and possible evolution effects in GRB relations (Li 2007; 
Tsutsui et al. 2008). However, it is found that no sign of evolution with redshift of the Amati relation, 
and the instrumental selection effects do not dominate for GRB relations (Ghirlanda et al. 2008, 2009). 
Nevertheless, further examinations of possible evolution effects and selection bias should be required 
for considering GRBs as standard candles for cosmological use. 
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